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ABSTRACT 

We observed IE 1207.4-5209, a neutron star in the center of the supernova 
remnant PKS 1209-51/52, with the ACIS detector aboard the Chandra X-ray 
observatory and detected two absorption features in the source spectrum. The 
features are centered near 0.7 keV and 1.4 keV, their equivalent widths are about 
0.1 keV. We discuss various possible interpretations of the absorption features and 
exclude some of them. A likely interpretation is that the features are associated 
with atomic transitions of once-ionized helium in the neutron star atmosphere 
with a strong magnetic field. The first clear detection of absorption features in 
the spectrum of an isolated neutron star provides an opportunity to measure 
the mass-to-radius ratio and constrain the equation of state of the superdense 
matter. 

Subject headings: pulsars: individual (IE 1207.4-5209) - - stars: neutron - 
supernovae: individual (PKS 1209-51/52) — X-rays: stars 



1. Introduction 

Although neutron stars (NSs) have been studied extensively for more than three decades, 
the properties of the superdense matter in their interiors still remain an enigma. We know 
neither the density nor the composition of a NS core. We are not even sure that NSs are 
indeed composed of neutrons — for instance, their cores could be pion or kaon condensates 
or a quark-gluon plasma. A key observational property that would help understand the true 
nature of these objects is the mass-radius relation — if we knew the masses and radii for 
a sample of NSs, we could compare them with the predictions of models based on various 
equations of state of the superdense matter, which are quite different for NSs of different 
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composition (see Lattimer & Prakash 2001 for details). A useful constraint can be obtained 
via measuring the gravitational redshift z of lines in the spectrum of thermal radiation 
emitted from the NS surface (atmosphere), which directly gives the mass-to-radius ratio: 
M/R = (c 2 /2G)[l — (1 + z)~ 2 ]. However, many attempts to detect spectral lines in thermal 
radiation from isolated (non-accreting) NSs have been unsuccessful. For example, no spectral 
lines have been found in the spectra of the Vela pulsar (Pavlov et al. 2001), anomalous X-ray 
pulsar 4U 0142+61 (Juett et al. 2002), and nearby radio-quiet NS RX J1865-3754 (Burwitz 
et al. 2001; Drake et al. 2002), despite sensitive observations with the Chandra grating 
spectrometers. In this Letter we report the first firm detection of absorption features in the 
spectrum of an isolated neutron star, IE 1207.4-5209. 

IE 1207.4-5209, a radio-quiet central source of the SNR PKS 1209-51/52 (also known 
as G296. 5+10.0), was discovered by Helfand & Becker (1984) with the Einstein observatory. 
Mereghetti, Bignami & Caraveo (1996) and Vasisht et al. (1997) interpreted the ROSAT and 
ASCA spectra of IE 1207.4-5209 as blackbody (BB) emission of T ~ 3 MK from an area 
with radius R ~ 1.5 {d/2 kpc) km. Zavlin, Pavlov & Trumper (1998) interpreted the observed 
spectra as emitted from a light-element (hydrogen or helium) atmosphere. For a NS of mass 
1.4 M Q and radius 10 km, they obtained a NS surface temperature T e g- = (1.4-1.9) MK and 
a distance d = 1.6-3.3 kpc, consistent with the distance to the SNR, d = 2.1+JJ kpc from 
the neutral hydrogen absorption measurements (Giacani et al. 2000). 

Zavlin et al. (2000) observed IE 1207.4-5209 with the Chandra X-ray Observatory and 
discovered a period of about 424 ms, which proved that the source is a NS. Second Chandra 
observation provided an estimate of the period derivative, P ~ (0.7-3) x 10~ 14 s s _1 (Pavlov 
et al. 2002a). This estimate implies that the characteristic age of the NS, r c ~ 200-1600 kyr, 
is much larger than the 3-20 kyr age of the SNR (Roger et al. 1998), while the conventional 
magnetic field, B e 3.2 x 10 19 (PP) 1/2 G = (2-4) x 10 12 G, is typical for a radio pulsar. 
Spectral analysis of these Chandra observations, which resulted in the discovery of two 
absorption features in the NS spectrum, is presented below. 

2. Observation and Data Analysis 

Chandra observed IE 1207.4-5209 on 2000 January 6-7 and 2002 January 5-6 with the 
spectroscopic array of the Advanced CCD Imaging Spectrometer (ACIS) in the Continuous 
Clocking (CC) mode. This mode provides the highest time resolution of 2.85 ms available 
with ACIS by means of sacrificing spatial resolution in one dimension. The source was imaged 
on the back-illuminated chip ACIS-S3. The effective exposure times of the two observations 
were 29.3 ks and 31.6 ks, detector temperatures —HOC and — 120 C, respectively. In both 
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observations the source-plus-background spectrum was extracted from a segment of 4" length 
centered at the source position; the background was taken from 20" segments in the one- 
dimensional images. The source count rates after background subtraction are 0.76 ± 0.01 s _1 
and 0.64 ± 0.01 s _1 in the first and second observations, respectively. The background 
contributes only about 3% of total counts, in the 0.4-5 keV energy band. The lower count 
rate in the second observation can be attributed to the reduction of the ACIS sensitivity at 
low energies with time, apparently caused by increasing deposition of a contaminant onto 
the ACIS filter 3 . 

The continuum models (e.g. BB, power-law, fully ionized NS atmosphere) failed to fit 
the source spectrum, which deviates very significantly from any of the models in a 0.5-2.0 
keV range, showing two absorption features near 0.7 keV and 1.4 keV. Another absorption 
feature observed near 2 keV remains uncertain due to the calibration uncertainty near the 
Ir M line from the telescope mirror coating, so we do not consider this feature further. We 
verified that the features at 0.7 and 1.4 keV are not due to uncertainties in the responses 
and background subtraction and concluded that they are intrinsic to the source. We also 
checked that no spectral features at these energies are seen in other Chandra observations 
with the same observational setup (e. g., observations of PSR B1055-52 and B0656+14). 

To characterize the absorption features, we use the spectral range from 0.4 to 2.5 keV 
and model the continuum as an absorbed BB with a temperature of kT = 0.26 keV. Since 
we do not know the exact nature of the observed features, we use three phenomenological 
models to check how sensitive the feature parameters are to the choice of model: 
lines produced by an absorbing layer with Gaussian profiles of absorption coefficients: 



/•'(/•:) /•; (/•:.) JJ ex P < ox P 

i=l,2 

multiplicative absorption lines with Gaussian profiles: 
F{E) = F C {E) J] {l-nexp 

i=l,2 ^ 



2r? 



(i) 



(E - Ejf 

2r? 



(2) 



lines with Gaussian profiles subtracted from the continuum: 
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In these models, F C (E) is the continuum spectrum, Ei and Tj are the central energy and 
Gaussian width of the i-th line, Tj and are the parameters characterizing the line depth (r^ 
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is the optical depth at the line center, and n — 1 — F(Ei)/F c (Ei) is the relative line depth, 
when the lines do not overlap, \E 2 — E ± \ >> Tj). 

An example of the fit to the January 2000 data with model (1) is demonstrated in the 
upper panel of Figure 1. The lower panel shows the continuum model from the upper panel 
and the residuals to demonstrate the profiles of the two absorption features. Including the 
lines in the spectral model improves the quality of the fit substantially: xt = 1-05 for 232 
degrees of freedom [d.o.f.] versus xl — 2.5 for 238 d.o.f. We find similar improvement in the 
fit quality for each of the line models for both observations. 

The fitting parameters for models (l)-(3) are presented in Table 1. The continuum 
normalization, R = (1.64 ± 0.02)(d/2.1 kpc) km, is essentially the same for all of the models 
for both observations. On the contrary, the hydrogen column density riu is substantially 
larger in the second observation. We interpret this as the result of an increased contamination 
of ACIS in the second observation which reduces the effective area at lower energies 3 and 
can be crudely modeled as an additional interstellar absorption (i.e., the % parameter does 
not represent the true interstellar absorption). 

The line parameters for models (1) and (2) are very close to each other within one 
observation. Model (3) results in somewhat lower central energies. The equivalent widths of 
the lines, calculated as W = J^[l- F(E) / F C (E)} dE (E min = 0.5 and 1.1 keV, £ max = 1.1 
and 2.5 keV, for the 0.7 keV and 1.4 keV lines, respectively), are virtually the same for the 
three models in a given observation. 

For a given model, the parameters of the 1.4 keV line are almost the same (within 
1.2er) for the two observations. The 0.7 keV line shows an apparent decrease of the central 
energy, width, and equivalent width. We attribute this to the effect of CCD contamination 3 
that results in an artificially large nu and distorts the parameters of the line. Since the 
contamination was presumably insignificant in January 2000, we believe that the parameters 
inferred from the first observation give more adequate description of the lines. 

For the empirical line models used, the best-fit models show residuals suggestive of 
asymmetric and/or multiple lines. Since the spectral resolution of the detector, about 100 
eV at these energies, is comparable to the model line widths Tj, each of the observed features 
can be fit equally well with multiple line models, with smaller intrinsic line widths. The 
observed broad features could also be produced due to a nonuniform magnetic field on the 
surface of the rotating NS. 

Reanalysis of the previous X-ray observations of this source has shown no features in 
the ROSAT/PSPC and ASCA/GIS spectra, which can be readily explained by the very poor 
spectral resolution of these proportional counters. The spectrum acquired with the ASCA 
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SIS (a CCD imaging spectrometer) does show features at the same energies, albeit at low 
significance due to insufficient statistics. The spectral models with two lines, which fit the 
Chandra spectra, are consistent with the ROSAT and ASCA spectra. 

3. Discussion 

Attempts to explain the absorption features as caused by the intervening interstellar 
(or circumstellar) material lead to huge overabundance for some elements. For instance, the 
absorption feature at 1.4 keV might be interpreted as a magnesium photo-ionization edge at 
E = 1.305 keV. The fit of this feature with a photo-ionization edge model gives the best-fit 
edge energy of 1.29 keV, not strongly different from the Mg K edge. However, to explain 
the strength of absorption, this model requires the Mg abundance about 600 times the solar 
value. The energy of the other feature, 0.7 keV, is close to the F K a (0.67 keV) and Mn L Q 
(0.64 keV) energies, but even larger overabundance of these elements would be needed to 
explain the observed feature. Therefore, the absorption in the ISM looks improbable, and 
the observed features are most likely intrinsic to the NS. 

There are two potential ways for these absorption features to be generated in the NS 
atmosphere — cyclotron lines and atomic transition lines. 

3.1. Cyclotron lines 

First we consider that the observed features are cyclotron lines produced in a strongly 
ionized NS atmosphere. If one assumes these are electron cyclotron lines, the features 
could be interpreted as the fundamental and the first harmonic of the electron cyclotron 
energy E ce = 1.16S n keV in a magnetic field B n = _B/(10 11 G) ~ 0.6(1 + z) or as two 
fundamentals emitted from two regions with different magnetic fields, B u = 0.6 (1 + z) 
and 1.2(1 + z), broadened by the radiative transfer effects and/or nonuniformity of the 
magnetic field. However, it is difficult to reconcile this interpretation with the expected 
strength of the surface magnetic field. The measured P and P imply a magnetic moment 
H ~ (3ic 2 PP/87r 2 ) 1/2 ~ 3 x lO 30 /^ 2 G cm 3 , where I = 10 45 / 45 g cm 2 is the moment of 
inertia. This magnetic moment corresponds to a magnetic field B e ~ 3 x 10 12 G at the 
magnetic equator, if the field is a centered dipole and the NS radius is R = 10 km. If the 
dipole is off-centered, or there are substantial multipole components, the maximum values of 
the surface magnetic field would be even higher. Although the inferred value of B e is based 
on P estimated from two observations under the assumption of a uniform slowdown (no 
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glitches, no timing noise, the NS is not in a binary orbit), the discrepancy of the fields looks 
uncomfortably large. In addition, the oscillator strength of the first harmonic is smaller than 
that of the fundamental by a factor of ~ E ce /(m e c 2 ) ~ 2 x 10 -3 (at E ce ^> kT — Pavlov, 
Shibanov, & Yakovlev 1980), so that it is hard to explain why the 1.4 keV feature is as strong 
as the 0.7 keV feature if we assume the two lines are associated with the same magnetic field. 

The magnetic fields of ~ 10 11 G should exist within the NS magnetosphere, and the 
resonant cyclotron scattering in such magnetic fields could absorb (scatter) ~ 1 keV pho- 
tons emitted from the NS surface (e.g., Rajagopal & Romani 1997). In a dipole magnetic 
field, B(r) — B e (1 + 3cos 2 #b) 1//2 (R/r) 3 , the resonant scattering of photons with energy E 
occurs at a distance from the NS surface r « 3.7R(B e /3 x 10 12 G) 1/3 (£/0.7keV)~ 1/3 (1 + 
3 cos 2 0b) 1 / 6 , where 9 B is the magnetic colatitude. To produce a line with a width 8E, 
the thickness of the scattering layer should be 8r ~ (1/3)(8E/E) r, i.e., 8r < 0.03r ~ 1 
km for the observed 8E/E < 0.1. It is not clear how two narrow layers (an analog of 
the Van Allen radiation belts in the Earth's magnetosphere?) could be formed. More- 
over, the electron number density needed to obtain a sufficient optical thickness r, n e ~ 
10 13 (r/0.5) (R/Wkmy 1 (£/0.7keV) 4/3 (B e /3 x 10 12 G)~ 1/3 (1 + 3 cos 2 6b)- 1 / 6 cm" 3 , signif- 
icantly exceeds the Goldreich- Julian density, n G j ~5x 10 cm" 3 , often adopted 
sonable estimate for the charge density in the NS magnetosphere. Therefore, this hypothesis 
does not look very plausible. 

Alternatively, one can assume that the spectral features are associated with ion cy- 
clotron energies, E ci = 0.63(Z/A)Bi 4 keV, where Z and A are the ion's charge and mass 
numbers. The surface magnetic field needed for this interpretation is > 10 14 G, much larger 
than the conventional magnetic field inferred for a centered dipole from the P, P measure- 
ments. However, if the magnetic dipole is off-centered to about 3 km below the surface, the 
magnetic moment // — 3 x 10 30 G cm 3 gives a surface field of about 10 14 G. We cannot ex- 
clude a possibility that the true period derivative is significantly larger than that estimated 
from the two observations (e.g., due to strong glitches), which would increase the inferred 
magnetic field. However, the ion-cyclotron interpretation of the observed features is not 
straightforward even in such a strong magnetic field. An explanation of two features as the 
fundamental and first harmonic of an ion cyclotron line is unrealistic because the ratio of 
the harmonic and fundamental oscillator strengths, ~ E ci /miC 2 , is extremely low. One could 
explain the two lines as produced by different ions, with (Z 1 A 2 ) / (Z 2 A±) = 2. For instance, 
one might assume that the 1.4 keV and 0.7 keV lines are due to protons (Z/A = 1) and 
alpha-particles (Z/A = 1/2), respectively, in a magnetic field B 14 = 2.2 (1 + z). The problem 
with this interpretation is that the strengths of the lines require similar abundances of hy- 
drogen and helium, which is hard to sustain because of element sedimentation in the strong 
gravitational field of the NS. To avoid this problem, one could assume that the radiation 
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emerges from a purely helium, partially ionized atmosphere — then the 1.4 keV line would 
be due to helium nuclei in B u = 4.4 (1 + z), while the 0.7 keV line would be caused by ion 
cyclotron transitions of once-ionized helium ions (Z/A = 1/4). The required field, however, 
is so strong that it is hard to expect a significant fraction of fully ionized helium at the effec- 
tive temperature of about 3 MK inferred from the spectral continuum (the temperature may 
be even lower if one applies NS atmosphere models). In addition, such cyclotron transitions 
should be accompanied by atomic transitions in, at least, once-ionized helium. 

3.2. Atomic lines 

The other possibility is that the observed features are atomic lines formed in the NS 
atmosphere. The available NS atmosphere models without magnetic field do not fit the 
observed spectrum, which is not surprising because the observed pulsations require a strong 
magnetic field. In recent years there has been significant work on structure and spectra 
of atoms in strong magnetic fields, mostly in fields below 10 13 G (e.g., Ruder et al. 1994; 
Mori & Hailey 2002, and references therein). Based on these works, we can exclude some 
possibilities. For instance, the observed absorption features cannot be explained as emerging 
from a hydrogen atmosphere because, at any magnetic field and any reasonable gravitational 
redshift, there is no pair of strong hydrogen spectral lines whose energies would match the 
observed ones. Therefore, one has to invoke heavier elements. A possible interpretation of the 
observed features as due to atomic transitions of once-ionized helium in a strong magnetic 
field, B ks 1.5 x 10 14 G, has been suggested by Pavlov et al. (2002b); detailed modeling 
of the spectrum will be presented in a forthcoming paper. This interpretation yields the 
gravitational redshift z = 0.12-0.23, which corresponds to R/M = 8.8-14.2 km Mq 1 . In 
that interpretation the 0.7 keV feature is a single line while the 1.4 keV feature may include 
several lines which cannot be resolved with ACIS. To confirm this, one should investigate 
the detailed structure of the detected strong features with a high-resolution spectrometer 
and detect other (weaker) helium lines. 

4. Conclusion 

The discovery of the absorption features in the spectrum of IE 1207.4-5209 provides the 
first opportunity to measure the mass-to-radius ratio and the magnetic field of an isolated 
NS. Measuring M/R is particularly important because it can constrain the equation of state 
of the superdense matter in the NS interiors, infer the internal composition of NSs, and test 
the theories of nuclear interactions. Our analysis of the low-resolution spectra has shown 
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that interpreting the observed features as cyclotron lines requires artificial assumptions, 
and such features cannot be explained as formed in a hydrogen atmosphere. The energies 
of the features suggest that they could be associated with the atomic transitions of once- 
ionized helium in an atmosphere with a strong magnetic field. This interpretation yields 
a gravitational redshift of about 0.17. To firmly identify the features and measure the 
gravitational redshift with better accuracy, deep observations with high spectral resolution 
are needed, such that would be able to resolve potentially multiple lines blended together 
due to the low spectral resolution of the CCD detector. 

We are grateful to Gordon Garmire, Leisa Townsley and George Chartas for the useful 
advice on the analysis of ACIS data. This work was partly supported by SAO grant G02- 
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NASA's Goddard Space Flight Center. 
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Table 1. Fitting parameters for three models for two observations of IE 1207.4-5209. 

nn,20 Ei ti Ti W\ £2 T2 T2 W2 

eV eV eV eV eV eV 

January 2000 

1 2.33 ± 0.44 750 ± 10 0.39 ± 0.03 125 ±13 108 ± 12 1432 ± 10 0.45 ± 0.04 105 ±12 102 ±14 

2 2.28 ±0.45 749 ±10 0.41 ± 0.03 134 ±14 112 ± 14 1432 ± 10 0.46 ± 0.04 112 ±14 104 ±13 

3 2.34 ±0.44 727 ± 12 0.40 ± 0.06 127 ±22 108 ± 15 1401 ± 12 0.43 ± 0.06 112 ±24 103 ± 13 

January 2002 

1 7.55 ± 0.54 735 ± 14 0.29 ± 0.03 94 ±20 64 ± 11 1423 ± 10 0.49 ± 0.04 125 ±14 129 ±14 

2 7.49 ±0.54 734 ± 14 0.30 ±0.03 100 ±21 68 ± 13 1422 ± 10 0.50 ±0.04 134 ±16 130 ±16 

3 7.48 ±0.54 721 ± 17 0.30 ±0.05 98 ± 36 66 ± 14 1380 ± 13 0.47 ±0.07 131 ± 26 129 ± 18 



Note. - - The "optical depth" n characterizes the line flux with respect to the continuum at the line center; 
for the second and third models it is defined as exp(— r^) = (see eqs. [2] and [3]). The errors quoted are the 
formal lcr uncertainties of the fits. The reduced xt values are within the range 1.03-1.06 (232 d.o.f. for the first 
observation and 236 d.o.f. for the second observation). 
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Photon Energy (keV) 
Fig. 1. — The observed count spectrum and the best fit with model (1) for the January 2000 
observation (top panel). The bottom panel shows the data and the continuum model. 



